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\ Chiral logarithms in m 2 are calculated at one loop, taking into account the 
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leading contributions to flavor symmetry breaking due to staggered fermions. 
I treat both the full QCD case (2+1 light dynamical flavors) and the quenched 
case; finite volume corrections are included. My starting point is the effective 



£> \ chiral Lagrangian introduced by Lee and Sharpe. It is necessary to understand 



the one- loop diagrams in the "quark flow" picture in order to adjust the 

calculation to correspond to the desired number of dynamical quarks. 
q ■ PACS numbers: 12.39.Fe, 12.38. Gc, 11.30.Rd 

I. INTRODUCTION 

> 

^ . Staggered (Kogut-Susskind, KS) fermions provide a competitive method for simulating 

Sh : 

full QCD, including the effects of virtual u, d and s pairs. Indeed, recent improvements in 
the action [l[ have made it possible to compute many physical quantities with rather small 
scaling violations in such "2+1 flavor" simulations |||3|]. 

The violation of KS-flavor symmetry, while reduced by the improved ( "Asqtad" ) action, 
remains however quite large at the lattice spacing (a m 0.13 fm) where most of the simula- 
tions in Refs. |2|U have been performed. The maximum splitting in mass squared among the 
various flavor pions, A max , is w (400 MeV) 2 at this lattice spacing. Because one can choose 
the staggered flavor on the valence quark lines, the flavor violations often enter through loop 
effects alone, and as such they have a typical size of only ~8% in most quantities. Here, I 
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estimate the size from a "typical" chiral logarithm, including flavor violation: 

^ln(A max /A 2 ), (1) 

where A is the chiral scale (taken, for instance equal to m p m 770 MeV), and f n = 131 MeV. 
Since the leading flavor violating terms in the improved KS action are 0(a 2 ), these dis- 
cretization effects can be reduced still more by going closer to the continuum limit; MILC 
simulations are currently in progress at a m 0.09 fm. 

On the other hand, if one focuses directly on chiral logs and works at a « 0.13 fm, 
the effect of KS-flavor violation should be large. Indeed, in Ref. 0, we were unable to fit 
m 2 /rri£ to the standard continuum chiral-log form M for 2+1 flavor QCD. (Here and below I 
stands for a generic light [u or d) quark; I neglect isospin violations throughout.) To test the 
hypothesis that KS-flavor violations are responsible for this behavior of m 2 /m£, one needs 
to compute the chiral logs in the presence of KS symmetry breaking. That computation is 
the subject of this paper f|. 

The effective chiral Lagrangian that describes KS fermions through 0(m, a 2 ) (m is a 
generic quark mass) has been constructed by Lee and Sharpe || . Their approach is to make 
a joint expansion in m and a 2 , which are considered to be comparably small. That is in fact 
the case for a simulation like 0, where the splittings in mass squared of various KS-flavor 
pions are comparable to the squared masses themselves. A one-loop calculations using the 
Lagrangian of Ref. |J would give the chiral logs (non-analytic terms) of 0(m 2 ,ma 2 ,a A ). 

The Lee-Sharpe Lagrangian corresponds to a single lattice KS field, which becomes 4 
degenerate flavors in the continuum limit. The main subtleties in the current work arise 
because I am interested in the more phenomenologically relevant 2+1 flavor theory being 
simulated, with m s ^ mp in general, not the 4-flavor (continuum degenerate) KS theory. I 
follow a three-step procedure: 

1. Generalize (almost trivially) the Lee-Sharpe Lagrangian to the 4+4 case, where one 
has two lattice KS fields with different masses (mp, m s ). This is an SU(8)l x SU(8)r 
theory, with mass terms and 0(a 2 ) KS-flavor violating terms breaking the symmetry. 
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2. Compute m% at one loop in the 4+4 case. I treat only the case where the pion is 
the U(1)a Goldstone particle. (This is the situation most immediately relevant to the 
fits attempted in [[J.) The symmetry then implies that m\ vanishes when — > 0, 
allowing one to simplify the calculation. 

3. Adjust the 4+4 answer to correspond to the 2+1 case of interest. This requires iden- 
tifying the contributions that correspond to £ or s virtual quark loops and multiplying 
them by 1/2 or 1/4 respectively. 

The subtleties are almost entirely in step 3. To identify those terms that should be 
reduced by a factor of 2 or 4, one needs to follow the "quark flow" approach, which was 
introduced by Sharpe to compute quenched chiral logarithms. Unfortunately, in the 
presence of KS-flavor symmetry violation, there does not appear to be an alternative, La- 
grangian approach ||[| that would automatically cancel the effects of the unwanted KS 
flavors. This corresponds to the fact that, in numerical simulations, the reduction from 4+4 
to 2+1 flavors is accomplished by taking the square root and fourth root of the I and s 
determinants, respectively. The procedure reduces each of the four KS flavors in a virtual 
quark loop equally, by a factor of 2 (4); it does not than cancel specific flavors. There is 
thus no equivalent (ultra) local lattice Lagrangian, and it is not at all clear how or whether 
one can represent the low energy chiral properties of this theory with an effective chiral 
Lagrangian [|l(| . 



Of course, in the continuum limit the KS-flavor symmetry is restored, and one may 
cancel any 2 of the 4 light quark flavors (and any 3 of the 4 strange quark flavors) with 2 
(3) bosonic pseudo-quarks. The effective chiral theory is then obvious: an "8|5" partially 
quenched chiral perturbation theory. This is discussed in Ref. |9| . Such a theory is trivially 
equivalent to standard 3 flavor chiral perturbation theory in the physical sector where the 
sea quarks are also the ones on the external lines. 

The remainder of this paper is organized as follows: Section |TJ presents the Lee-Sharpe 
Lagrangian and generalizes it to the 4+4 case. The calculation of to one loop is then 



described in Sec. |T|. In Sec. [TV], the quark flow picture is applied to compute in the 
2+1 case. The analytic terms at 0(m 2 ,ma 2 ,a 4: ), necessary to perform consistent fits that 
include the chiral logs, are also touched upon. The quenched case follows in Sec. [V]. Finite 
volume corrections are then discussed briefly in Sec. [VT]. Finally, I make some remarks about 
the extension of the current approach to other physical quantities in Sec. |V11[ That section 
also contains some comments about the fits of MILC data to the chiral forms derived here; 
however, a detailed discussion of such fits and their consequences will be presented elsewhere 

n . 



II. LEE-SHARPE LAGRANGIAN FOR 4+4 FLAVORS 

Lee and Sharpe |J first construct the continuum effective Lagrangian for KS lattice 
fermions, including all terms at 0(a 2 ). The KS-flavor symmetry is broken by 6-quark 
operators. They then write down the corresponding effective chiral theory to leading order 
in a 2 and m, the quark mass. 

For a single KS field (4 continuum flavors) the Lee-Sharpe chiral Lagrangian describing 
the pseudo-Goldstone bosons has a (nonlinear) S't/(4)^ x S77(4)# symmetry, broken by the 
mass term and by the KS-flavor violating operators. The 4x4 matrix E is defined by 

16 

E(x) = exp(i<f>/f), ct> = Y,(t>aT a (2) 

o=l 

where cj) a are real, / is the tree-level pion decay constant (I use the normalization /„. = 
131 MeV), and the Hermitian generators T a are chosen as follows, 

T a = {£5, i^B, i^u, 1} ■ (3) 

Here ^ are the flavor gamma matrices, £^5 = £^£5, ^ = (with fi < v in eq. @), and 
I is the 4x4 identity matrix. S transforms by S — > under SU{4)l x SU(A)r. 

Note that I include the singlet meson (pi, corresponding to generator J, in eqs. (0) and 
(0). The mass of the ?/-like meson (often called $o> U P to overall normalization) gets a large 
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contribution (= mo) from the anomaly in the Lee-Sharpe case and is not included in their 
formalism. But I want to generalize to two KS fields (8 continuum flavors). In that case, 
only the SU(8) flavor singlet is heavy and can be integrated out. 577(4) singlets, e.g., a 
singlet "kaon" made of one quark from each KS multiplet, can be light. I therefore leave 
the singlet as well as the m term explicit in intermediate steps, but will ultimately take 
m -> oo.Q 

The (Euclidean) Lee-Sharpe Lagrangian is then 

f 2 1 ?m 2 

£ (4) = ^tr(^S^St) - -pnftrp + £t) + -^{hf + a 2 V , (4) 

where \i is a constant with units of mass, the coefficient of (fii is conventional,^ and V is the 
lowest order KS-flavor breaking potential: 

- V = d tr(&£&£ + ) 

+C 2 i[tr(S 2 )-tr(e 5 Se 5 S) + /i.c.] 
+C7 3 ±]T[tr(£,£^£) + / i . C .] 

V 

+c4J2M^^) + h. c .} 

V 

+C 5 \ ]T [tr(e.Se.Et) - tr(U^B^)\ 

V 

+C e X; [tr(^S^Et) . (5) 

fl<V 



1 Sharpe and Shoresh [12 show that integrating out the <3?o is mathematically equivalent to keeping 
it in the calculations and taking the mass parameter mo — > oo at the end. If that limit is to be 
taken, it is unnecessary to include $o dependence other than the niQ term in the action. However, 
in the quenched case (Sec. |V|), other <3?o dependence will need to be considered. 

2 mo/6 by definition multiplies ((fin + 022 + ^33 + • • ~) 2 ■ Here, the normalized (fix = 5 Ya=i fin- 
This definition of m 2 , is independent of the number of flavors and corresponds to that in Refs. HHJ 
(after renaming /i 2 — ► ttiq) and that in Ref. @ (after choosing N = 3 there). With the current 



definition of /, the parameter 5 introduced in Q is equal to m^/ {24it 2 f 



2\ 
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Effects of order a 2 m or a 4 are neglected in eq. (|J). 

The potential V has an "accidental" 5*0(4) flavor symmetry (evidenced by the sums 
over v and ft), which is larger than the lattice symmetry group. The pions fall into 5 5*0(4) 
representations with flavors £5, £^5, £ M , I. Expanding eq. (f|) to quadratic order in 0, 
the pion masses are found to be: 



AiTfl 2 

ml B = 2/im + -^Sb,i + a 2 A^ B ) , (6) 

where, in a convenient abuse of notation that will be used from now on, the index B takes 
the 16 values {5, fib, \w{p, < v),fi,I} and £7 = /. The m\ contribution to the mass of (pi is 
shown explicitly, and A(£b) is given by 

A(£ 5 ) = 

A(£ M5 ) = ^(Ci + 2 + 30 3 + 4 - O s + 30 6 ) 
A(^) = ^(20 3 + 20 4 + 40 6 ) 



A(U = ^(Ci + 2 + 3 + 30 4 + O s + 30 6 ) 



A(/) = ^(40 3 + 40 4 ) (7) 

A(£ 5 ) vanishes because of the staggered-flavor nonsinglet Ua(1), which is represented in the 
chiral theory by 

£ -> e* & £e* & , (8) 

with 6 1 the £7a(1) angle. Since this symmetry is unbroken by the lattice regulator, 7r 5 is a 
true Goldstone boson in the chiral limit. 

From simulations such as those in Ref. 0, one learns that the 4 term in eq. (j5|) is 
the largest contributor to the 0(a 2 ) flavor violation ||. This leads to approximately equal 
splitting between the pions, in the order 7r 5 , 7t m5 , tt^, n^, The other operators are not 



3 Omitting the disconnected terms in the 717 propagator (as is often done in simulations - see, 
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entirely absent, however; their contributions to the splittings are of order 10% of that of C4 
fTT| . The reason for the smallness of Cj, i 7^ 4, is not known.Q 

We now need to generalize to the case of two KS fields with different masses me and m s , 
i.e., 4+4 continuum flavors. The field E in eq. (0) becomes an 8 x 8 matrix, given by 

/ 7T K\ 

E(x) = exp(i$//) $ EE , (9) 

\tf s) 

where the 4x4 fields tt, K, and S describe "pions," "kaons," and "ss" mesons, respectively. 
As in eqns. (HIH), n = S a =i ^a^oj with 7r a real, and similarly for K and S, except that the 
K a are complex. The SU(8) singlet oc tr7r + trS will be eliminated below by taking the 
m — > 00 limit. 

Defining the 8x8 mass matrix A4, 

/m/J \ 

M =[ , (10) 

V m s l) 

the SU(8)l x SU(8)r Lagrangian that generalizes £(4) of eq. @ is then 

f 2 1 9m 2 

£ (4+ 4) = V tr (^d M £t) - 7 /i/ 2 tr(A4E + MSt) + — ° (tt/ + 5,) 2 + a 2 V . (11) 
o 4 o 

Generalization of the KS-flavor breaking potential V in eq. ([5]) requires a little thought. 
For a single KS field, the symmetry breaking 4-quark operators in the effective continuum 
theory have the generic form 

Q(7s®60QQ(7s'® ( 12 ) 
where 75 is an arbitrary spin matrix, £p is an arbitrary flavor matrix, and the indices 
S, S', F, F' are contracted in various ways determined by the lattice symmetries. The effec- 
tive chiral operators may then be found from the 4-quark operators by treating £p and £p' 
as spurion fields 0. 



e.g., j3j), eliminates the mo contribution to the singlet mass and extends the approximately equal 
splitting rule to KS flavor /. 
4 More precisely, the approximately equal splittings imply only that Cq, C\ + C2, and 2C3 — C5 
are small. I thank M. Golterman for pointing this out. 



With two KS fields on the lattice, there is an exact vector SU(2) (broken only by mi ^ 
m s ) that mixes them|] The symmetry guarantees that the KS-flavor breaking 4-quark 
operators now have either of the two forms: 

Qids ® &)Qj Qj(is> ® &>)Qi , Qiiis ® Z F )Qi Qj(is> ® 6«)Qj , (13) 

where i and j (= 1,2) are SU(2) indices. The operators in eq. ( |I"3l) are "flavor mixed" 
and "flavor unmixed," respectively. By Fierz transformation, the flavor mixed operators 
can be brought to the flavor unmixed form, so we may assume all 4-quark operators are 
of the latter type.0 But if all 4-quark operators are flavor unmixed, then all spurion fields 
are SU(2) singlets. In other words, for V in eq. (|ll|) we may take simply eq. (|5]) with the 
replacement 

(Zb \ 

6* - • (14) 

VO 

The tree level masses of the pions, kaons, and S (ss) mesons are then 



m\ B = 2fim e + a 2 A(£ B ) 
m K B = ^i. m i + m s) + a 2 A(£e) 

m% B = 2fim s + a 2 A(£ B ) , (15) 

where A(£b) are given by eq. (0), and the m terms are not included because we will treat 
2m\{j[i + Si) 2 /3 in eq. (PD as a vertex (summed to all orders) below. 



5 There are also corresponding axial symmetries, whose generators are the direct product of the 
L^4(l) generator, 75 (g) ^5, with the vector SU(2) generators, but I will not need them here. 

6 This may be at the expense of mixing the color indices (suppressed in eqs. (12) and (|T3|)), but 
they are irrelevant in the corresponding chiral effective theory because of confinement. 
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III. ONE-LOOP PION MASS FOR 4+4 DYNAMICAL FLAVORS 



The calculation of m 2 to one loop in the 4+4 theory (eq. ([H])) is now straightforward. 
I confine my attention to corrections to the Goldstone pion mass m WB . The cases of the 
non-Goldstone pions and the kaons, which will also be useful in confronting simulation data, 
will be left to future publications. 

The graphs are all tadpoles, with vertices coming from the kinetic energy term 
tr((9^S(9 M E^), the mass term tr(.ME + .ME^), and the symmetry violating term V. They are 
shown in Figs. [II ^ and |^, respectively. 

I write the one-loop 7r 5 self energy as[| 

S (^ 2 ) = 9^7^ Sl + S 2° n + S ^ C ) , (16) 

where I have explicitly separated the contribution coming a disconnected propagator in the 
internal loop, E^ (corresponding to Fig. 0(b)), from the other, connected contributions in 
Ej 011 . The terms "connected" and "disconnected" are applicable at the QCD level: the dis- 
connected diagrams contain gluon intermediate states in the meson loop (Fig. f|(d),(e),(f)). 
The Ua(X) symmetry guarantees that the sum of E™ 11 and Ejj is proportional to m 2 B . 
However, this is not true for individual diagrams contributing to Ef 311 , so I do not include 
the m 2 factor in the definition. 

For the moment, I calculate only the universal chiral logarithms, and ignore the analytic 
terms in the diagrams (which means the divergences are also ignored) . The logarithms come 
from the following integrals: 

Ji s / wji -> i^ m2 lnm2 (17) 

_ f d 4 q q 2 2 f d 4 q 1 4 2 

1 2 = / - — — — = — m 1\ + / - — — — > -m mm . (18) 

J 2vr 4 q 2 + m 2 J 2vr 4 16n 2 V ; 



7 The distinction between the field matrix E in the Lagrangian and the self energy E - both 
conventional notations - should be clear from context. 
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Si gets contributions only from Fig. |](a). The result is 

Ei -> -8m^ 5 lnm^ 5 - 8m^ lnm^ - £)m^ fl lnm^ , (19) 

B 

where no sum on /x is implied, B takes the 16 values {5, /x5, < z/),/i, J} as usual, and 

the arrow means that only the chiral logarithms are included, as in eqs. ( |TTD and (|T8|). 
Figure |I|(b) contributes only to E| 011 : 

- 8m** In m^ 5 + 8m 4 ^ In + £ < a In ™L > ( 20 ) 

again with no sum on \x . Note that among the pions, only the vector and axial vector flavors 
contribute in eqs. (|I^) and (pop. It is easy to show, using the fact that [£5,!] = [£5,6/^] = 
[£5^5] — 0, that no terms of the form 7r 4 , tt^tt 2 , or n 2 ^] can be generated by the kinetic 
energy term. The absence of coupling to irj explains why the kinetic energy term does not 
contribute to E^ . 

Figure ||(a) contributes only to Ef 311 . Couplings of 775 to all the other pions are generated 
by the mass vertex, so there is no restriction on the terms that can enter. I find 

EnTjAqc — > — 3mi \ia ml — Ami mi him 2 

2, MADS 7T5 7T5 7T5 7T M 5 7T M 5 

-lSm^m^lnm^ -4m^m^lnm^ [no sum on /1, 1/] 

-3ml m 2 ^ ln ""4, ~ i m l 5 + m UE m L lnm L • ( 21 ) 

Here, the factors of mjj. or mj <5 in each term come from the quark masses, which are 
proportional to Goldstone particle masses squared by eq. (|15|) with A (£5) = (eq. (0)). 

The mass vertex also generates the only contributions to the disconnected diagram, 
Fig. 0(b). This is the sum of 1, 2, 3, ■ • ■ insertions of the vertex ^^(717 + Sj) 2 on the internal 
717 line. The 717 line in turn comes from a term proportional to tt 2 tt] in the mass vertex; 
the same term generates the connected contribution —3m\ i m 2 Ki ln(m^ ) in eq. (PT|). A single 
insertion of the m\ vertex, which would be the only contribution in the quenched case, gives 

v disc,quench _ Q 2 f ^1 4mo/3 
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The standard quenched double pole is evident. Iterating the Wq vertex and summing the 
geometric series then results in 

S disc = 3m 2 f A 4m 0/ 3 m) 

[q 2 + ml) 2 1 + -2T%r- + — 



There is no here need to take the m\ — > oo limit or to perform the integral; the conversion 
to the 2+1 flavor case can best be done directly from this expression. Note that 717 and Si 
intermediate states now enter. As discussed in Sec. |V|, these correspond to virtual £ and s 
quark loop contributions. 

Finally, we need to consider the KS-flavor violating vertices, eq. (|5|). It is not hard to 
show that, as in the kinetic energy case, only the vector and axial pions contribute. From 
Fig. [3], I obtain (no sum on //) 

E§$ - -8a 2 A(£, 5 )< 5 lnm; 5 - 8a 2 A(^)m^ lnm 2 M - £ a'AfoX* ln ™L , (24) 

B 

where A is given in eq. (0). 

We can now put the results together to find the 7r 5 mass at one loop. Writing 

K -loo P) 2 = ^ + _i_ e5 (25) 

we have 

_ yicon 1 v^con , v^con , ^disc 2 /qj-n 

e 5 — Zj 2,KE "T Zj 2,MASS ^.V ^2 m, 7r 6 Zj l ' V Z0 / 

From eqs. (pj) , (|2~0|), (pl|), and fl24|) we then have 



e 5 -> ^' 5 (-3m 2 5 lnm 2 5 + I2m 2 ^ 5 lnm 2 ^ - 18m^ lnm 2 ^ 



+12m^ him^ - 3m^ lnm^j + [no sum on i/] . (27) 

where the arrow, as usual, means that only the chiral logarithms are kept. Note that e§ is pro- 
portional to m 2 5 , as it must be by Goldstone's theorem. (E^ already has a factor of m 2 5 in 
eq. (p3|).) The term J2b m K B m ( m /f s ) ( ec L- (POD ) combines with — J2b a2 ^(^s) m l' s m ( m i: s ) 
(eq. (0)) to cancel against —m 2 Kb ^ B m 2 KB \n(m 2 KB ) (eq. fl2T|)). Here I have used 
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m 2 KB =m 2 K5 +a 2 A^ B ) , (28) 
which follows from eqs. (|15|) and (0). Similarly, the terms 8m^ g ln(m^ 5 ) and 8m^. ln(m^) 



in eq. ([H]) combine with corresponding terms in eq. (|24T ) to produce terms that have a 
factor of m\ h . Indeed, eq. (p4|), the effect of the symmetry violating vertices, is completely 
determined by the requirement that the full answer go like m 2 5 (up to logs) in the chiral 
limit. This makes it clear, for example, that the Ci coefficients in eq. (|3p must enter eq. ( |24| ) 
only through the combinations A(£b), and that the final answer can depend only on the 
various meson masses. 

Terms of the form m\ h m\ B \\\{w? KB ) also cancel between eqs. (|19D and (|2l|) . This is 
"accidental" in the sense that it is not required by symmetry. It corresponds to the fact 
that there happen to be no \n(m 2 K ) terms in the standard continuum 2+1 flavor result [|J. 
For similar quantities, such as f n or (OImmIO), the \rt(m 2 K ) terms will not cancel 0]. 

Finally, note that in the symmetry limit (m WB = m n B = m %fiy = m n ^ = m ni ) the entire 
one-loop correction in eq. (^) is proportional to E^ . This is also true of the standard 
result H; the reason will be explained in the next section. 



IV. ONE-LOOP PION MASS FOR 2+1 DYNAMICAL FLAVORS 

To adjust the result in eq. (|27D to the case of 2+1 dynamical flavors, we need to identify 
the contribution from each of the possible quark flow diagrams shown in Fig. ^. Many of 
the arguments used to identify the quark flows will be familiar from Ref. 0. Once the 
contributions from virtual quark loops are determined, the adjustment to 2+1 dynamical 
flavors is accomplished by multiplying every £ (i.e., u,d) quark loop by 1/2 and every s 
quark loop by 1/4. 

First note that Figs. f§(d), (e), and (f) have disconnected internal meson propagators, 
and correspond to Fig. 0(b). Figure ^(d) has a single ml insertion and therefore generates 
the contribution in eq. (P2"|). Figures |^(e) and (f) represents the iteration of the m\ vertex 
through the introduction of virtual quark loops. Together, Figs. [|(e) and (f) change eq. (|2~2D 
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to eq. (H). 

To identify contributions from the remaining graphs in Fig. [|, we first study the pos- 
sible 2 into 2 meson scattering diagrams in the quark flow language, Fig. [5[ Letting 
i,j,k,n = 1,2, ••■8 be quark flavor indices, Fig. ||(a) corresponds to a meson vertex of 
the form $jj$jj$fc ra $ n fc, where $ is the meson field in eq. Similarly Fig. ||(b) corre- 
sponds to $jj$jfc$fc n $ n i. Now, consider the case where all four quark flavors i,j,k,n are 
different. Then vertices generated by the kinetic energy or mass terms in eq. (|TT|) cannot 
have the form $jj$jj$fc n $ n jt because they are formed only by a single trace of fields (and 
the diagonal matrix M). Since Fig. |5|(a) cannot depend on the flavors of the quark lines to 
this order (except trivially through A4), this implies that Fig. |(a) is not generated at tree 
level by the kinetic or mass terms for any values of i, j, k, n. 

The KS-flavor violating vertices behave differently. Because at least some of the £b 
explicitly entering eq. (|5]) must be off-diagonal, it is not hard to see that Fig. ||(a) can in 
fact be generated. Note, however, that since are diagonal under the vector 577(2) that 
mixes the two KS fields, Fig. |5|(a) still vanishes when i,j are of £ type (upper 4x4 block) 
and k, n are of s type (lower 4x4 block). 

Now consider Fig. |4](c). This graph is absent when the quarks in the virtual loops are 
of s type because the needed vertex vanishes to this order. If we then choose m s = me, it 
can make no difference whether a quark in a virtual loop is £ or s type, so this graph must 
also vanish when the loop quarks are of £ type. A slight generalization of this argument 
shows that the graph again vanishes when there is one virtual quark of each type (i.e., a 
kaon loop).0 

Once we know that Fig. ||(c) is absent, intermediate kaons could only be generated by 

8 Add a third (s') KS quark field, plus a pseudoquark KS field s' of the same mass. This partially 
quenched theory is then equivalent to the original theory in the I, s sector. The graph Fig. |||(c) 
then vanishes when there is one s and one s' in the loop. Now choose m s / = ni£. 
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Fig. §(b). Since the K contribution to 65 in fact vanishes to this order, graphs of type 
Fig. 0(b) with an s type quark in the virtual loop must cancel among themselves. Once 
more choosing m s = me, this also implies that Fig. 0(b) cancels for a i type virtual quark 
loop. 

The surviving connected contribution to €5 comes from Fig. 0(a). Furthermore, since 
the vertex in Fig. |](a) is of the type in Fig. |^(a), which can only be generated by KS-flavor 
violating terms, such connected contributions must vanish in the symmetry limit. This is 
in fact the case, as was mentioned at the end of the last section. Because Fig. 0(a) contains 
valence quark lines only, these contributions do not change when we go from 4+4 to 2+1 
dynamical flavors. 

The only required adjustment is therefore to the disconnected contribution, S^ 1550 . We 
need to divide every virtual quark loop contribution in Figs. |(e) and (f) by 2 or 4, de- 
pending on whether it is t or s type, respectively. In eq. (|23|) this is easily accomplished 
by the replacements (4ml/3)/(q 2 + m 2 7 ) — ► (2ml/3)/(q 2 + m 2 j ) and (Aml/3)/(q 2 + m 2 Sl ) — ► 
(mo/3) /(q 2 + rrtg ) in the denominator. In other words, we have in the 2+1 case: 

E *= = 3m 2 (tl 4m o/ 3 (29) 



We can now take the limit m\ — > 00 and put eq. fl2Tj| ) in the form 

Sf C = <f^[ {q2 + mli) - {q 2 + m 2j) > ( 30 ) 

where m 2 = (2m| 7 + m 2 i )/3. Equation ([[7]) then gives 

£f c -> m 2 5 (6m 2 , In m*, - 2m 2 qi In m 2 m ) . (31) 

Inserting this into eq. (p7|), we get the 2+1 result for the chiral logarithms 

e 5 — > m 2 (-3m 2 mm 2 + 12m 2 In ml r - 18m 2 lnm 2 

+12m 2 ti In m 2 ^ + 3m 2 1 In m\ i — 2m 2 7 lnm 2 \ [no sum on ji, u] . (32) 



In the KS-symmetry limit, this gives the standard result [0. As discussed above, the result 
comes entirely from E^ 1 * 30 m this limit. 
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What about the analytic terms? These can in general come from a quite complicated 
set of operators in the chiral theory: standard 0(m 2 ) ("p 4 ") operators 0(ma 2 ) (e.g., 
Vtr(.MS + .ME^)) or 0(a 4 ) (e.g., V 2 ) chiral operators generated by the same terms in the 
continuum effective action that led to eq. flTT|), and entirely new chiral operators of 0(ma 2 ) 
and 0(a 4 ) coming for example from terms in the continuum action that Lee and Sharpe |3| 
call Sq F( " B \ which break the symmetries of eq. (|TT|) down to the lattice symmetries. However, 
as a function of quark mass, the analytic terms in (m*~ loop ) 2 of the relevant order can only 
be proportional to m,, m 2 or m,m s , since they must vanish as m, — > 0. This means that 
terms of 0(a A ) in the joint m, a 2 expansion cannot enter. 

Putting eq. fl32|) in eq. (|25|), adding in the analytic terms just discussed, and writing the 
overall m 2 5 as 2/xm, (eqs. (|T5|) and ([FD), I arrive at the final result in the 2+1 case (no sum 
on ji, u): 

{1 / m 2 m 2 m 2 

1 + 16^( - 3< ^ + 2< 5 In ^ - Sm 2 ^ In 

+2m 2 M In + |m 2 z In _ ln + 2me K 3 + (2m, + m s )K' 4 + a 2 C\, (33) 

where A is the chiral scale and K 3 , K' A , and C are independent of m and a to this order. The 
term a 2 C may alternatively be considered as a discretization correction to \i. K' A is related 
to the parameter of Ref. Q by 

ir 4 = (2m, + m s )if; . (34) 

In the KS-symmetry limit eq. (|33[) reduces to the result in ||. 

Since the term a 2 C is what absorbs mass- independent cutoff dependence in eq. (j33"D, 
we can make a rough estimate of its size by computing the change of the logarithms when 
A changes between two reasonable values, say 0.5 and 1.0 GeV. For MILC simulations at 
a « 0.13 fm, for which A max w (400 MeV) 2 (see discussion before eq. ([]])), I find a 2 C ~ 0.05. 
Discretization corrections to K 3 and K f 4 come from higher, 0(m 2 a 2 ), contributions not 
considered in eq. (|3"3"D. I expect that such corrections will be significant because they can 
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be generated by the same operators in the effective QCD Lagrangian that produced the 
KS-symmetry breaking potential, V. 



V. ONE-LOOP PION MASS IN THE QUENCHED CASE 

Given the understanding of the quark flows in Fig. |] developed in the previous section, the 
result in the quenched case is easy to write down. E!^ in eq. (|27|) changes to the quenched 
version, v^ lsc > c i ucnch ) w hich corresponds to Fig. |](d) only. The other terms in eq. (|27|) , which 
arise from Fig. §(a), are unaffected. We thus have, in the quenched case 



e 5 — > mi. (-3m 2 lnm^ + YlmL , \arni , - 18mi In 



mi 



+12m 2 ^ him 2 ^ — 3m\ i lam 2 A + £jdisc,quench ^ Q gum Qn ^ ^ ^ 

•j-U vidisc, quench • 1 /F7ol\ 

with L 2 given by eq. (|22|). 

Since wig is not to be taken to infinity in the quenched case, it is necessary to consider 
other <pj dependence in the Lee-Sharpe Lagrangian,! 9 ] eq. (Q). In Ref. || it was shown that 
it is sufficient to add in a contribution to the </>/ kinetic energy: (2a/3)(d fl (f)i) 2 . This merely 
changes ml in eq. ( f22f ) to m\ + aq 2 . 

The non-analytic terms in vjd i&c >q ucnch ma y Iqq extracted using 

f d A q 1 -13 1,2 

2 3 = - — — — — — — = — - — Ji — > -mm (36) 

3 7 (2vr) 4 (q 2 + m 2 ) 2 2m 9m 16tt 2 V ; 

J4 " / (2^ (^bp = 2^^ J2 = Jl " m2j3 leb 2m2lnm2 ' (37) 

where Xi and X2 are defined in eqs. (|T7D and (p~8|). 

Putting the logarithms in vjd isc -q uench together with eq. (|35|), and using eq. (p5|), gives the 
final quenched answer: 



9 In the quenched case, there is no point in going to the 4+4 Lagrangian, eq. (11), since s type 
quarks do not couple to the pion. 
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{1 / 771 ^ TTl^ TV? 

1 + 16^ ( - ln # + 2 < ln - 3m '- ln 

+2< In ^ - 1< ln ^ - ln ^ + §o< In ^) + 2m^ 3 + a*c} , (38) 

where of course the parameters fi, K% and C may be different from those in eq. fl3~3"|). I 
keep the term multiplying a separate from the other lnm^ 7 term, since they arise from 
different integrals and therefore will have different finite volume corrections. Analytic terms 
depending on or a are not included since they can be absorbed in /x and A. Note that 
the parameter K' 4 does not appear here because it multiplies m s , which cannot enter in the 
quenched pion mass. 

In the symmetric limit, eq. (]38[) agrees with results in || or 0. As usual, the power 
counting in the quenched case is somewhat problematical; one must assume that Wq and a 
are in some sense small to stop at one loop ||[7||. 



VI. FINITE VOLUME CORRECTIONS 

The technique for computing the finite volume corrections to the results of Sees. [TV] and 
|V] is standard |13| , ^4| . However, I will provide enough detail here to allow the reader to 
include the finite volume corrections numerically in fits to eqs. ([33|) and (^). I assume 
that the Euclidean time dimension, T, is large enough that it may be taken infinite, and 
only corrections from the finite spatial dimensions (length L) need be considered. The 
generalization to the case where T is also finite will be obvious. 

The results in finite volume come simply from replacing the integrals Zj (i — 1, . . . , 4) in 
eqs. flT7j), (|T8|), ( p6|) and (|37|) by their finite volume counterparts T\ L \ We have, for example, 



= Iv/ d ^ 1 (391 

L^J (2v)ql + {q n y + m? 1 ' 



q n = 2-nn/L ; n = (n x , n y , n z ) (40) 

where n x , n y , and n z are integers, and I am assuming periodic boundary conditions. 
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4 L) and zj L) are denned analogously. 

Since the integrals of interest are divergent, it is preferable to work with the differ- 
ences Z$ — Xj, which are finite as long as care is taken in introducing a cutoff. Examples 
of acceptable regulators include multiplying all integrands/summands by (A 2 /(g 2 + A 2 )) fc 
(where A; is a large enough power to render the integrals convergent) fT5|l , point-splitting the 
corresponding position-space propagators, and defining the theory on a lattice. 

Assuming the integrals/sums have been regulated appropriately, I define the dimension- 
less finite volume corrections Si by 

(41) 
(42) 
(43) 
(44) 





-Xi = 


m 




[mL) 


16tt 2 


T (L) 


-x 2 = 


4 

m 


8 2 


[mL) 


167T 2 


T (L) 


-x 3 = 


1 




[mL) 


16vr 2 


T (L) 


-x 4 = 


2 

m 


Si 


[mL) 


16vr 2 



The correction 8\ can easily be written in terms of the Euclidean position space propa- 
gator, defined in infinite volume by 

where K\ is a Bessel function of imaginary argument. The corresponding finite volume 
propagator, is given by 

G<1> <*> = h ? / ^y gww = G <*> + £ G <* + M > < 46 » 

with q = (q , q n ). The last equality follows from the Poisson resummation formula, or simply 
by noting that it is the unique solution of the corresponding differential equation with the 
correct (periodic) boundary conditions. Subtracting eqs. ( ^6|) and (^5|), setting x = 0, and 
putting the result into eq. fl4"l|) using eqs. (^) and (|i~7|) gives 

(5i(mL) = — ^ K i(\n\ mL ) (47) 
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For practical values of ml, 8\ can be readily computed to machine precision with this 
formula. (One is limited of course by the accuracy to which K\ is evaluated.) An alternative 
approach, inspired by Ref. |16|], is to compute the difference of the (regulated) sum and 
integral directly in momentum space, treating the sum as an integral by rounding Lq/(27r) 
to the nearest n and using a standard numerical integration program |T7[] . With a reasonable 



regulator (e.g., the one from |15 |, mentioned above), this method gives results consistent 
with eq. (f47|), but it is much slower. 

The other finite volume corrections follow from the relations among the integrals given 



in eqs. (ITq) , (|3q) and (p7|) and the standard recursion relation K[(z) = —K\{z)jz — K (z), 



where prime implies differentiation. I find 

5 2 (mL) = -8x(mL) (48) 
5 3 (mL) = -8 x (mL) - —5[(mL) = 2 £ K (\n\mL) (49) 

2 n^O 

5 A {mL) = 25 1 {mL) + f ^—5' l {mL) = S^mL) - 5 3 (mL) (50) 
In deriving eq. ((48|), one must be careful to regulate so that the difference between the sum 



and the integral of unity vanishes in the limit of infinite regulator. The regulator in [15 
works well, for example. 

Large mL asymptotic expressions for the Si are obtained by including only the 6 nearest 
neighbors in the sums in eqs. ( flTD and ( |49"D and using the only leading terms in K\ and Kq. 
I find 

8 1 {mL) = -8 2 {mL) ~ ^?e~ m£ (51) 

(mL) 2 



5 3 (mL) ~ -5AmL) ~ Q x ^- e - mL . (52) 

V mL 

While these expressions are useful for estimating the size of finite size effects, the leading 
corrections to eqs. ( |51~D and ( |52|) have just one higher power of l/(mL). In the quantitative 
applications I have in mind ]T1|, for which mL is in the range 3 to 9, the full expressions 



(eqs. ([470 , (f48D, (fl9|) and (p0|)) — or at least the full K\ and i^o from the nearest neighbors 
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— will be needed. It is important to note, however, that for very small quark masses, mL<l, 
eqs. (|47|), (pE8D , ( f49|) and (^) are not applicable, no matter how many terms are included 
in the sums. In that limit, the zero-mode (the spatially independent component of the pion 
field) must be treated exactly |T4[ . 

Comparing eqs. pi|)-(pH|) with eqs. flT7|), (|T8D, ( |3T)| ) and ( |3"T| ) then gives the rules for 



including finite volume corrections. In the 2+1 case, one needs merely to make the replace- 
ment 

2 2 
Tfl Tfl 

In— ^ln— + 5 1 (mL) (53) 

for every logarithm in eq. (§3|). It is interesting to note that when rri£ — > (with L and 
the symmetry violating parameters, Cj, fixed), the leading finite size correction to eq. ( |33|) 
is negative, due to the — |m^ 5 lnm^ 5 term. In the standard case, the leading correction is 
positive, due to the positive sign of m\\&m\. 

In the quenched approximation, the replacements in eq. ( |3~8D are: 

ui 2 

tit? m? 

a ln — > a ln + \a,8i{m^jL) (54) 

2 2 
Tfl Tfl 

In — — - — > ln — — + 5\ (mL) [all other logarithms] . 
A 2 A 2 



VII. REMARKS AND CONCLUSIONS 

My main results are given for the 2+1 flavor case in eq. (|33|) and for the quenched 
case in eq. (j38|). Finite volume corrections to these formulas appear in eqs. (|53| ) and (|54"|), 
respectively, with the values of 5i given in eqs. ( fi?|) to (pTp. 

I have computed only 77i^ 5 , but a generalization of this calculation to vn? Kh is straight- 
forward and is in progress flIB|. Slightly more complicated, but also on the list of things 



to compute [113] are /„., fx, (^^), and the non-Goldstone 7r and K masses. These quanti- 



ties do not vanish in the chiral limit, so there will be fewer short-cuts available. A further 
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complication for a quantity like f w is that there is no accidental cancellation of diagrams 
corresponding to Fig. f|(b). An additional adjustment for the number of flavors in virtual 
quark loops will therefore be required. However, the adjustment is not difficult, since him 2 K 
contributions come uniquely from this diagram. To generate the 2+1 case, one will merely 
have to: (1) multiply the lnm K terms by 1/4 to count the s quarks in the loop correctly, 
and (2) subtract half the (original) lnm K terms, after putting m s — > me, to count £ quarks 
in the loop correctly. For quantities such as fx or m K , where there will be valence \nm K 
contributions, the In m\ and In m| terms (from the connected internal propagator) can be 
used to make the adjustment. 

Preliminary fits of eq. fl33|) to MILC 2+1 data are encouraging: the fits give good con- 
fidence levels for reasonable ranges of quark mass | 1T| . The contrast with the standard fits 
that do not take into account KS flavor breaking || is especially striking. I emphasize here 
that the new fits have the same number of free parameters as the standard ones. This is 
because the various pion splittings are not left free in eq. (|33| ) but are determined first by 
the data for the masses of all the pions. 

In the quenched case, one generally treats m (or equivalently b = m\j (247r 2 / 2 )) as a free 
parameter. As such, decent fits can apparently be performed whether or not one takes flavor 
breaking into account. However, the values of 5 obtained seem to be considerably closer to 
the real world values (5 ~ 0.18) when flavor breaking effects are included. In both the 2+1 



and quenched cases, we are collecting additional lattice data at small quark mass |TT| , which 
should significantly increase the discrimination between the symmetric and flavor-breaking 
fits. 
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V 



(a) 



(b) 



FIG. 1. Chiral perturbation theory graphs contributing to the pion propagator from kinetic 
energy vertex (solid square). The external lines are Goldstone pions, i.e., 7r$. The dots represent 
the derivatives in the vertex. In (a) they act on the external lines; in (b), the internal. 




V 




(a) 



(b) 



FIG. 2. Same as Fig. ||, but from the mass vertex (solid triangle). The internal propagator in (a) 
is the connected propagator only (no tbq insertions), even when it is neutral (717). All disconnected 
contributions are in (b); i.e., the cross represents one or more insertions of the m§ vertex. 
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FIG. 3. Same as Fig. [T], but from the flavor breaking vertex V (solid circle). 
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(a) (b) 




(e) (f) 

FIG. 4. 

Quark flow diagrams that could contribute to the (flavor non-singlet) pion propagator 
at one loop in the chiral expansion. Diagrams (d), (e) and (f) are disconnected graphs that 
correspond to figure Fig. 0(b); (f) represents the sum over two or more intermediate virtual 
quark loops. Note that (d), (e) and (f) each have a second version, in which the two sides of 
the meson loop come from different valence lines. For the purposes of this paper, one may 
consider the diagrams shown to be generic, subsuming both versions. Diagrams with "$ 
vertices," such as Fig. 3(i) of Ref. J7|, are not relevant: either because they are eliminated by 
the m\ — > oo limit (the 2+1 case), or because they have virtual quark loops (the quenched 
case). Diagrams with odd numbers of <3> lines at a vertex, such as Fig. 3(f) of Ref. can 
be eliminated by a field redefinition ||. 
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